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COPYRIGHT NOTICE 

[01] Portions of the disclosure recited in this specification contain material that 

is subject to copyright protection. Specifically, source code instructions by which 
specific embodiments of the present invention are practiced in a computer system are 
included. The copyright owner has no objection to the facsimile reproduction of the 
specification as filed in the Patent and Trademark Office. Otherwise all copyright rights 
are reserved. 



Cross- References To Related Applications 

[02] This application is related to the following co-pending U.S. Patent 

Applications which are hereby incorporated by reference as if set forth in full in this 
specification: 

Serial No. [TBD], filed on [TBD], entitled "SINGLE INSTRUCTION 
MULTIPLE DATA IMPLEMENTATIONS OF FINITE IMPULSE RESPONSE 
FILTERS"; and 

Serial No. 10/057,694, filed on January 23, 2002, entitled "METHODS 
FOR EFFICIENT FILTERING OF DIGITAL SIGNALS." 
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Background Of The Inv ntion 

Field of the Invention 

[03] This invention is related in general to computer processing and more specifically 
to the use of single instruction multiple data (SIMD) instructions to achieve finite 
impulse response filter operations in a digital processor. 

Description Of The Background Art 

[04] Finite Impulse Response (FIR) filter operations are an important type of digital 
computation or processing. FIR filters are commonly used, for example, in pre- 
processing, post-processing, motion compensation, and motion estimation for video 
compression standards. The implementation of FIR filters in computer programs, or 
other digital processing approaches, is useful in many other applications including audio 
processing, signal conditioning, simulation of electronic components, etc. 
[05] FIR filter operations can be very demanding on digital processing systems 
because of the large number of iterative operations that must be performed very quickly. 
The number of operations, speed of operation, resolution of coefficient values, and other 
factors all contribute to the accuracy of the implementation and the amount of processing 
resources that are necessary to achieve a design goal. In this respect, a slight advantage 
in FER filter operations that are executed frequently (i.e., in an "inner loop" of a program) 
can result in very significant performance gains. 

[06] An FIR filter that is of special interest in video compression and encoding 
techniques is referred to as a transversal or tapped delay filter. These filters multiply a 
set of coefficients to pixel values of a video frame to generate a new pixel value. Such an 
operation is useful, for example, to compress an image by combining adjacent pixel 
values into a smaller number of pixel values. Typically, this type of FIR filter includes 
only positive coefficients. 

[07] Fig. 1 illustrates four pixel values a u a 2 , a^ and a A . Subpixel b is desired to be 

the average of the four pixels computed as: 

[08] b = (a x +a 2 +a 3 +a 4 +2)»2, (1) 
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where » is a bitwise right shift operator. 

[09] In a typical application where pixel values are limited to values in the range 0- 
255, the pixels a lv . .,a 4 are each represented in one byte or 8 bits. Thus, a total of four 
bytes is necessary to operate on the four pixel values at once. 

[10] Typically, a single frame in a digital video presentation of moderate resolution 
can include 600x800 = 480,000 pixels. Such a frame might be displayed 30 times per 
second. Moreover, it may be necessary to perform additional "passes" over the frame so 
that, for example, in subsequent passes the pixels, themselves, are combined into 
subpixels to further compress an image. Thus, numerous subpixel computations may be 
necessary. Further digital video formats, such as high-definition television, use much 
higher screen resolutions and color depths. It should be apparent that such filter 
operations could place enormous requirements on processing resources, especially when 
the operations must be performed in real time. 

[11] One approach that the prior art uses to provide increased efficiency in filter or 
array operations is to use Single Instruction Multiple Data (SIMD) instructions. Such 
instructions allow value-packing, byte-packing, or other concatenating of values into a 
single word or other unit of data. The unit of data can be processed quickly by 
performing a desired operation in parallel on the packed values. 

[12] SMD-type instructions are available in many processors. Examples include Intel 
Multi-Media Extensions (MMX)™ and Streaming SIMD Extension (SSE)™, as well as 
NEC VR5432, Equator MAP-CA™, and Philips TM-1300 processors. In processors 
whose architecture supports SIMD instructions there are typically multiple identical 
processors, AT, each with its own local memory where it can store data. All processors 
work under the control of a single instruction stream issued by a central control unit. 
There are typically N data streams, one per processor. The processors operate 
synchronously: at each step, all processors execute the same instruction on a different 
data element. This architecture allows N computations in parallel. Thus, if N=8, it is 
possible to achieve a computational speedup of 8. 

[13] Fig. 2 provides an example of the operation of a SIMD instruction. The SIMD 
instruction performs an operation, "OP," on two sets of data: A=[a u .. .,a 8 ], a vector of 8 

data values, each of which is an unsigned 8-bit integer, i.e., ^€[0,255]; and 5=[6 b . . .,6 8 ], 
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another vector of unsigned integers within the range [0,255]. The final result C=[c 1? . . .,c 8 ] 
is achieved by simultaneously operating on all 8 values of a { and b i as c t = a { OP b i9 for 
i=l ,. . .,8. In this example, A, B and C are 64-bit registers in which all 8 values of a i9 b i9 
and c { are packed as contiguous bytes as shown in Figure 2, i.e., JV=8. Such operations are 
also known as packed operations, since 8 values of data are packed in a single register^, 
BorC. 

[14] One specific type of operation of interest in filter operations is the PA VG 

operation that can be found, e.g., in the Intel MMX™ instruction set. The PAVG 

instruction performs the following computation: 

[15] PAVG(A,B) = [(ai+bi+1) » 1, z=l,...,8]. (2) 

[16] This operation takes 8-bit values of a,-, b h and stores the intermediate sum 

1) in 9 bits before doing bitwise logical right shift operation to get the final result. 
It is available in many processors, including the ones mentioned above, and uses only one 
instruction. This instruction has the latency of 1 clock cycle in the Intel Pentium III, 2 
clock cycles in Intel Pentium 4, and Advanced Micro Device's (AMD's) Athlon, with a 
throughput of 1 clock cycle. The same performance is realized for other operations in 
these architectures, such as packed addition (+), subtraction (-), bitwise AND (&), bitwise 
OR (|), bitwise EXCLUSIVE-OR ( A ), bitwise right shift (»), and bitwise left shift («) 
operations. 

[17] Although SIMD instructions can improve the efficiency and speed of 
computations, such instructions are sometimes difficult to use effectively when the SIMD 
instructions do not provide the exact type of operation needed. For example, as stated 
above, PAVG computes (a^b^l) » 1. An average of two vectors rounded up . 
However, it is more desirable in some filter operations to obtain which is a 

truncated average where the remainder, or fractional part, is discarded. Such a difference 
in operation is significant where multiple passes of frame data are made as the average 
intensity value of subpixels may increase and result in artifacts or other objectionable 
qualities to the processed data. In the architectures discussed herein, a SIMD instruction 
to compute (<*,•+&,•) » 1 is not provided. Typically, a non-SIMD approach must be used. 



4 



[18] A problem also arises when the number of arguments required by a SIMD 
operation is not the same as the number of variables in a formula to be implemented by 
the SIMD operation. For example, if a SIMD instruction accepts two arguments then it is 
"mismatched" to implement a formula, computation or operation with more than two 
variables or values. The same can be said, for example, for a SIMD instruction with 
three arguments used to implement a formula with other than three variables, etc. 
[19] Fig. 3 illustrates a non-SIMD approach to compute (a i +6 l )»L 
[20] In Fig. 3, a i9 b t are unsigned integers within the range [0,255], i.e., each a t , b t is 
represented in 8 bits. The number of processors, N=8, i.e., the operation (a,-+ft/)»l is 
simultaneously performed on 8 values of a i and b i for i— 1,. . .,8. All 8 values of a i 
(usually contiguous pixels) are packed in 64-bit register, A, and 8 values of b t in 64-bit 
register B. Since a i +b i can exceed 8 bits, the 8-bit (byte) values of a i9 b t are unpacked into 
16-bits (words) as four 16-bit values per 64-bit register,. Then the packed registers A and 
B are added together, followed by bitwise logical right shift by 1, followed by packing 
again. Note that in most processors, data can be packed into 64-bit registers as 8 (byte), 
16 (word), 32 (dword), or 64 (qword) bit values only. Figure 3 shows the conventional 
method of doing the packed operation c l -=(a / +6 / )»l for z=l,. . .,8. It is clear from Figure 
3, that given sufficient memory, 9 instructions are needed to achieve the result 
Cf=(fl/+fc/)»l for all 8 values of a t and Each instruction in Figure 3 is represented by 
an ellipse. 

Summary of Embodiments of the Invention 

[21] The invention provides improved results in some cases of digital calculation of 
finite impulse response (FIR) filters. A preferred embodiment of the invention is applied 
to techniques for FIR calculation discussed in the co-pending patent application entitled 
"SINGLE INSTRUCTION MULTIPLE DATA IMPLEMENTATIONS OF FINITE 
IMPULSE RESPONSE FILTERS " referenced, above. In the co-pending patent 
application a system for efficient derivation of FIR values is presented using single- 
instruction multiple data (SIMD) types of operations. In a preferred embodiment, the 
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results of the FIR calculations are subjected to additional operations using a SIMD 
instruction called PA VG. The results of PA VG are a rounded-up average of two sets of 
packed values. Adjustments are made on the rounded-up average to obtain an exact 
desired result for various filter calculations, or to obtain results within a desired error 
range, or results that do not exceed, or fall below, desired values in relation to the exact 
answer. Various techniques for minimizing processor resources (e.g., processing cycles, 
memory) are presented. 

[22] These provisions together with the various ancillary provisions and features which 
will become apparent to those artisans possessing skill in the art as the following 
description proceeds are attained by devices, assemblies, systems and methods of 
embodiments of the present invention, various embodiments thereof being shown with 
reference to the accompanying drawings, by way of example only, wherein: 
[23] One embodiment of the invention provides [@@] 
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Brief Description of the Drawings 

[24] Fig. 1 illustrates a subpixel average of four pixel values; 

[25] Fig. 2 shows an example of the execution of a single-instruction multiple-data 

(SIMD) instruction; 

[26] Fig. 3 shows a non-SIMD approach to a calculation; and 
[27] Fig. 4 shows a SIMD implementation with adjustment. 

Detailed Description of Embodiments of the Invention 

[28] The technique of the present invention includes adjusting an FIR calculation 
result using SIMD instructions to obtain an improved result. This technique is the focus 
of section 5 of this specification. Other sections include text from the co-pending patent 
application entitled "SINGLE INSTRUCTION MULTIPLE DATA 
IMPLEMENTATIONS OF FINITE IMPULSE RESPONSE FILTERS," cited, above, 
upon whose results the adjustments of the present invention are based. 
[29] A preferred embodiment of the invention uses Intel's MMX/SSE architecture, 
including the SIMD PA VG operation. Other embodiments may use other processors, 
instructions and operations in a manner similar to that disclosed herein and realize similar 
computational benefits. In addition, other techniques and approaches for performing 
processing may benefit from one or more of the features presented herein, such as the 
techniques of the related patent application "METHODS FOR EFFICIENT FILTERING 
OF DIGITAL SIGNALS," cited above. 
[30] Table I shows notations used in this application. 



Operator 


Description 


+ 


Addition 




subtraction 


& 


bitwise AND 


1 


bitwise OR 


A 


bitwise exclusive OR 
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» 


bitwise logical ri*?ht shift 


« 


bitwise logical left shift 




Bitwise NOT 


CLIP(x) 


Clips x to range [0,255] 


ODD(x) 


Returns 1 when x is odd, 0 otherwise 


EVEN(x) 


Returns 1 when x is even, 0 otherwise 



TABLE I 



[31] The present invention allows computing c l =(fl l +6 / )»l for i=l,. . .,8, in an efficient 
manner using a SMD instruction such as PA VG. Note that simply using the PA VG 
instruction on packed values in registers A and B will not yield the correct answer. For 
example, when afrb t is an odd number PA VG(a i9 b t ) gives a result that is one more than 
the correct answer. The result of a PA VG operation must be adjusted as follows: 

C = PAVG(A, B) - (A A B) & 0x01, 
where 0x01 is a 8-bit number whose least significant bit is 1 and the rest are 0's. 
[32] The PA VG operation with adjustment is shown in Fig. 4. Assuming sufficient 
memory, only 4 instructions instead of the previous 9 instructions (without using PAVG) 
are needed to achieve the packed operation C=(A+B)»l. This is an approximate speedup 
of9/4 = 2.25 times. 

[33] A preferred embodiment of the invention achieves the same computational result 
as in Fig. 4 with even fewer instructions by appropriately using the PA VG instruction in 
combination with supplemental logical operations to adjust for the rounded-up average. 
As described below, several FIR filtering operations can be modified to obtain result in 
fewer instructions when compared to conventional SMD implementations. 
[34] Without loss of generality, let A x , A 2 , . . . A \ 6 be 1 6 vectors, each of which contain 
8 packed data elements. For example, A 5 contains 8 data elements A 5 = • . a (5,8)l- 

Each data element a^ l6 i) for z=l,...8, is within the range [0,255], i.e., they are 

represented by bytes, and A h . . ty A ]6 are packed 64-bit registers: 

A r K/,i)> • a (/,8)] for/=l,...,16. (4) 
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[35] We perform various operations on the packed 64-bit registers A x , . . . , A x 6 to obtain 
different FIR filters described below. We define packed 64-bit vectors/registers ONE and 
ONE 4 as follows: 

ONE= [0x01, 0x01, 0x01, 0x01, 0x01, 0x01, 0x01, 0x01], 
ONE 4 = [0x0001, 0x0001, 0x0001, 0x0001], (5) 
where 0x01 is a byte containing 1 in its least significant bit and 0's elsewhere. The 
packed 64-bit register ONE contains 8 packed bytes, each containing 0x01 . On the other 
hand, the packed 64-bit register ONE 4 contains 4 packed words (16 bits), each containing 

0x0001. 

[36] The FIR filters used in a preferred embodiment include: 

1. Type 1 Filter: (A l +A 2 +c*ONE) » 1, where ce {-2,-1,0,1,2}, (6) 

2. Type 2 Filter: (A l +A 2 +A 2 +A 4 +c*ONE) » 2, where ce {0,1,2}, (7) 

3. Type 3 Filter: (A x +A 1 +A 2 +A A +A s +A 6 +A 1 +A%+c*ONE) » 3, where ce {0,1,2,3,4}, (8) 

4. Type 4 Filter: (A x +A 2 +...+A X5 +A X6 +c*ONE)» 4, where c€ {0,1,2,3,4,5,6,7,8}. (9) 

[37] All 4 types of FIR filters are useful for video compression applications. There are 
numerous FIR filters that can be constructed from these 4 basic types, in addition to those 
described herein. For example, the filter (2A x +A 2 +A 3 +2*ONE) » 2 is a Type 2 filter 
with ^=^4. Similarly, the filter (A x +2A 2 +2A 3 +2A 4 +A 5 +4*ONE) » 3 is a Type 3 filter 
with A 2 = 1 Afr A 7> =A 1 , and A A -A % . Many other types of filters can be constructed as will be 
apparent to one of skill in the art. 

[38] Instructions according to the present invention can be used to obtain exact filter 
computations. Such exactness may be necessary as, e.g., in motion compensation and 
estimation applications where accuracy is key. In other cases an approximation of the 
filter computation may be sufficient. For example, in cases where the number of 
operations is large an approximate computation can be a better tradeoff. The 
approximations of the preferred embodiments produce an error of ±1 in the final result 

for a small percentage of all values of a^e [0,255] for i=l 8, and 7=1, ...,16. These 

results are useful in cases such as post processing, where a small error of ±1 (in intensity 
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or color value) is inconsequential in the final result. Naturally, other approximations of 
different degrees of accuracy are possible and are within the scope of the invention. 

L Type 1 FIR Filters 

[39] There are 5 variations of the Type 1 FIR filters (A } +A 2 +c*ONE), where 
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ce {-2, -1,0,1,2}, based on the 5 choices of constant c. We state the SIMD 

implementation for each of these filters: 

(A x +A 2 -2*ONE) » 1 =PAVG(A l9 A 2 ) - ONE- (A X A A 2 ) & ONE, (10) 

(A x +A 2 -ONE) » 1 = CLIP(PAVG(A X A 2 ) - ONE), (11) 

(A x +A 2 )» 1 = PAVG{A x A 2 )-{A x a A 2 )8lONE, (12) 

(A x +A 2 +ONE)» 1 = PAVG(A X ^ 2 ), (13) 

(A x +A 2 +2*ONE) » 1 = PAVG(A X A 2 ) + K^z) & QAffi). (14) 

[40] There is a less efficient solution for (A x +A 2 ) » 1 that will be used to simplify 

expressions: 

(A x +A 2 ) » 1 = (A x » 1) + (A 2 » 1) + (A x &A 2 & ONE). (15) 
[41] Although (15) uses more instructions that (12), we need this expression to 
evaluate other filters. In (15), (A x & A 2 & ONE) is a correction term that is necessary 
when both A x , and A 2 contain odd integers. An approximate solution for (A x + A 2 ) » 1 is: 
(A x +A 2 ) » 1 = PAVG(CLIP(A x -ONE)4 2 ) or PA VG(A x ,CLIP(A 2 -ONE)). (16) 
[42] In most processors, subtract and CLIP() can be realized in one instruction. So the 
implementations in (16) require only 2 instructions. 

II. Type 2 FIR Filters 

[43] There are 3 variations of Type 2 filters (7) based on the 3 choices of constant c, 
where ce {0,1,2}. We show the derivation of each filter. We define the following 64-bit 
packed registers, each containing 8 data elements of one byte each: 

B x =PAVG{A X A 2 \B 2 = PAVG{A Z A A \EB X ={A^ (17) 

A. Type 2, Filter 1 : R = (A i + A-> + A? + A 4 + 2*ONE) » 2 
i. Conventional SIMD Solution 

[44] This filter can be implemented in SIMD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (4 Instructions) A iL = Unpack Low 4 Bytes of A h for i = { 1 ,2,3,4} , 

2. (4 Instructions) A m = Unpack High 4 Bytes of A b for i = { 1 ,2,3,4} , 
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3. (5 Instructions) Add and Shift lower 4 words of A x ,...^i 4 to obtain lower 4 words of 
i? L as: 

R L = (A 1L + A 2L + ^ 3L + ^ 4L + 2*CW£ 4 ) » 2, 

4. (5 Instructions) Add and Shift higher 4 words of A y ,...^i 4 to obtain higher 4 words of 
J? H as: 

R H = (A m + A 2U + A m + A 4H + 2*ONE A ) » 2, 

5. (1 Instruction) Pack R H and R L into final register R. 

We require 19 instructions to perform this filter by conventional SIMD methods. 
ii. Efficient SIMP Solution 

[45] In order to implement this filter efficiently, we simplify as follows: 

R = ({{A x +A 2 +ONE)»X) + ((A 3 +A 4 +ONE)»l) + E)»\= {B x +B 2 +E) » 1, (18) 
where E is the correction term that is necessary when both (A x +A 2 +ONE) and 
(Aj+Aj+ONE) are odd integers as in (15). Detection of odd or even integers is performed 
with the functions ODDQ and EVENQ. Where ODDQ returns "1" for each packed 
argument value only if the packed argument value is an odd number and returns "0" 
otherwise, and where EVENQ returns "1" for each packed argument value only if the 
packed argument value is an even number and returns "0" otherwise. 

E=ODD{A x +A 2 +ONE) & ODD(A 3 +A A +ONE) = EVEN(A X +A 2 ) & EVEN(A 3 +A A ) 

= ~(A^A 2 ) & ~04 3 M 4 ) & ONE = ~{EB X \ EB 2 ) & ONE. (19) 
We note that Ee {0,1}. From (18) and (12), we have: 

_{PAVG(B x ,B 2 )-(B l A B 2 )&ONE when£ = 0 
~{ PAVG(B U B 2 ) when£ = l* 

We simplify (20) as: 

R = PA VG(B X ,B 2 ) - (B { A B 2 ) &~E& ONE, 

which is same as: 

R = PA VG(B X ,B 2 ) - (B X A B 2 ) & ((A X A A 2 ) | (AfA 4 )) & ONE. (21) 
The solution in (21) requires 10 instructions. We have an approximate 19:10 (approx. 
2:1) speedup by using (21). 
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iii. Approximate SIMP Solution 

[46] Besides the accurate solution, we can obtain an approximate solution in fewer 
instructions by assuming the least significant bit of EB X or EB 2 as 0 or 1. Assuming the 
least significant bit of EB X or EB 2 = 1 , we get: 

R s PA VG(B X ,B 2 ) - (B X A B 2 ) & ONE. (22) 
This solution requires 6 instructions, and according to (16), it is close to the following: 

R = PA VG(CLIP(B x -ONE),B 2 ) or R = PA VG(B x ,CLIP(B 2 -ONE)) . (23) 

This solution requires only 4 instructions, and produces a maximum error of ±1 in the 
final result for 12.5% of all possible values of A h ...,A4 between [0,255]. The error never 

exceeds ±1. We get a computational efficiency of 19:4, nearly 5 times speedup. 

B. Type 2. Filter 2; R = (A t + A-> + A, + A, + ONE) » 2 
i. Efficient SIMP Solution 

[47] As seen in Section 3. A, this filter can be implemented by conventional SIMD 
methods in 19 instructions. For efficient implementation, we simplify as follows: 

R = (((Ai+A 2 +ONE)»\)+((A 3 +A 4 )»\)+E)»l = (B x +B 2 +(E-{EB 2 &.ONE))) » 1. 

(24) 

Here EB 2 is due to the correction term in (12), and E is the correction term in (15) as: 
E = ODD(A x +A 2 +ONE) & ODD(A 3 +A 4 ) = EVEN(A X +A 2 ) & ODD(A 3 +A 4 ) 

= ~(A , A A 2 ) & (A 3 A A 4 ) & ONE = (~EB l & EB 2 ) & ONE. (25) 
We note that E T = (E -(EB 2 &ONE)) e {0,-1}. From (24), (1 1), and (12) we obtain: 

_(PAVG(B u B 2 )-(B ] A B 2 )&ONE when£ r =0 
{ PAVG(B u B 2 )-ONE when£ r =-f 

Note that (E-(EB 2 &ONE)) = -1 when (A { A A 2 ) & (A 2 A A 4 ) & ONE = 1. We simplify (26) 
as: 

R = PAVG{B { ,B 2 ) - {{B X A B 2 ) \ {{A X A A 2 ) & (A 3 A A 4 ))) & ONE. (27) 
The solution in (27) requires 10 instructions, an approximate 19:10 (nearly 2 times) 
speedup. 
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ii. Approximate SIMP Solution 

[48] We can obtain four approximations of (27) by assuming the least significant bit of 
EB\ or EB2 as 0 or 1 . A good approximate solution is with the assumption that the least 
significant bit of ZTi?2=0, which gives us the same solutions as (22) and (23), which 
require 4 instructions and has a maximum error of ±1 for 12.5% of all possible values of 
A\,..., A A g [0,255]. We get a computational advantage of 19:4. 

C. Type 2s Filter 3: R = (A 1 + A^ + A, + A* ) » 2 

i. Efficient SIMP Solution 

[49] The filter can be implemented by conventional SIMD methods in 17 instructions. 
For efficient implementation, we simplify as follows: 

R = (((^i+^ 2 )»l)+((^ 3 +^ 4 )»l)+E)»l = (B l +B 2 +(E-(EB } +EB 2 )&ONE))»\. 

(28) 

Here EB X and EB 2 are due to the correction term in (12), and E is the correction term in 
(15) as: 

E=ODD(A { +A 2 ) & ODD(A 3 +A A ) = (A { A A 2 ) & (A 3 A A 4 ) & ONE = EB X & EB 2 & ONE. 

(29) 

We note that Ef^iE-iEB^EB^&ONE) e {0,-1}, and R is same as (26). Note that E T = - 
1 when (A X A A 2 ) | (A 3 A A 4 ) & ONE = 1. We simplify (26) as: 

R = PAVG(B l9 B 2 ) - ((*, A * 2 ) I {A X A A 2 ) \ (AfAJ) & ONE. (30) 
The solution in (30) requires 10 instructions. We have an approximate 17:10 speedup by 
using (30). 

ii. Approximate SIMD Solution 

[50] We can obtain four approximations of (30) by assuming the least significant bit of 
EB\ or EB 2 as 0 or 1. A good approximate solution is with the assumptions that the least 
significant bits of EB\ or EBi =1, which gives us: 

R = CLIP(PAVG(B h B 2 ) -ONE). (31) 
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This solution requires 4 instructions and has a maximum error of ±1 for 12.5% of all 
possible values of ^ lv ..^ 4 e[0,255]. We have a computational advantage of 17:4, approx. 
4 times. 

D. Type 2, Special Filter 1: R = (2A i + A^ + A i + 2*ONE) » 2 

[51] This filter is same as Filter 1 with A x = A 2 . It can be implemented by conventional 
SIMD methods in, e.g., 17 instructions. This type of filter is used extensively in, for 
example, standards proposed by the Joint Video Team (JVT) as, for example, in 
[CHANGE THIS - CHANCHAL TO UPDATE TO MORE CURRENT 
REFERENCE -> ISO/IEC MPEG and ITU-T VCEG, Geneva, Switzerland, Oct., 02; 
entitled "Editor's Proposed Draft Text Modifications for Joint Video Specification (ITU- 
T Rec. H.264 | ISO/IEC 14496-10 AVC), Geneva modifications, draft 26.and other 
coding schemes. <r END REFERENCE] 
[52] We can simplify (21) as: 

R=PAVG(A U B 2 ) - {AfB 2 ) & (A 3 A A 4 ) & ONE. (32) 
This solution requires 7 instructions. One can verify that (32) is close to the following: 

R = PA VG(A U PA VGiCLIP(A 3 -ONE)A 4 )), (33) 

which requires only 3 instructions instead of 17 instructions by conventional SIMD 
methods, a nearly 6 times speedup. However, (33) produces an error of ±1 for a very 
small 0.1% of all possible values of A u ... 9 A 4 €[0,255]. 

E. Type 2, Summary of Results 

[53] Table II below summarizes the instructions required to compute each filter (given 
sufficient memory) by the efficient and conventional SIMD methods. For the Efficient 
case, we give the instructions required for the exact and approximate solutions. 



Summary of Results for Type 2 FIR Filters. 



Type 2 Filters 


Conventional 


Efficient Method 


Speedup 


Method 


Exact 


Approx. 


Exact 


Approx. 
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(4i + A z + A 3 + A A + 2*0NE) » 2 


19 


10 


4; - • 


,1.9 


4.75 


(A l +A 2 +A J +A 4 + ONE)»2 


19 


10 ■ ■ 


• :: ;4>. .. 


1.9 


4.75" . 


{A x +^ 2 +^ 3 +^ 4 )»2 


17 


h: 10 




1.7 


4.25 


(Z4, + A 2 + A 3 + 2*ONE) » 2 


17 




3 


2.4' 


5.67 



TABLE II 

[54] The shaded areas show significant improvements in efficiency due to the analyses 
developed here. 



3. Type 3 FIR Filters 

[55] There are 5 different Type 3 FIR filters depending on the 5 choices of c in (8). We 
define the following packed 64-bit registers, each containing 8 data elements of one byte 
each: 

B x = PAVG(A l9 A 2 ) 9 B 2 = PAVG(A 3 4 4 ), B 3 =PAVG{A 5 A 6 \ B 4 =PAVG(A 79 A S ) 9 
C x = PAVG(B h B 2 ), C 2 =PAVG(B 3 >B 4 ) 9 
EB X = (AfA 2 ) 9 EB 2 = (A 3 *A 4 ) 9 EB 3 = (A 5 *A& EB 4 = (A 7 A A$) 9 

EC X = (B^B 2 ) 9 EC 2 = (BfB 4 ). (34) 

A. Type 3, Filter 1: R = (A i + A 2 + + A A + A g + A* + A „ + A 0 + 4*ONE) » 3 
i. Conventional SIMP Solution 

[56] This filter can be implemented in SMD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (8 Instructions) A iL = Unpack Low 4 Bytes of A i9 for i = 1 , . . . ,8, 

2. (8 Instructions) A m = Unpack High 4 Bytes of A i9 for / = 1,. . .,8, 

3. (9 Instructions) Add and Shift lower 4 words of ^ 1? ...^4 8 to obtain lower 4 words of 
# L as: 

R L = (A , L + A 2L + A 3L + A 4L + A 5L + A 6L + A 1L + A SL + 4*ONE 4 ) » 3, 

4. (9 Instructions) Add and Shift higher 4 words of A } v ..y4 8 to obtain higher 4 words of 
R H as: 
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R H = (A m + ^t 2H + ^ 3 H + ^4H + ^5H + ^6H + ^7H + ^8H + 4*OA® 4 ) » 3, 

5. (1 Instruction) Pack i? H and i? L into final register R. 

We require 35 instructions to compute this filter by conventional SIMD methods. 



il New SIMP Solution 

[57] In order to implement this filter without unpacking, we simplify it as follows: 
R = (((A } + A 2 + A?+ A 4 +2*ONE)»2) + 

((A s + A 6 + A 7 + ^ 8 +2*ONE)»2) + E) » 1 , (35) 
where E is the error/correction term that is necessary for dividing up the expression into 
two parts each containing a Type 2, Filter 1 (21). Note that according to (21) and the 
registers in (34), we have: 

(A x +A 2 +A z +A A +2*ONE)>>2 = C X - (E x & ONE), 
(A 5 +A 6 +A 7 +A s +2*ONE)»2 = C 2 - (E 2 & ONE), (36) 
where E x = EC X & {EB X \ EB 2 ) and E 2 = EC 2 & (£5 3 \ EB 4 ) are error/correction terms 
obtained in (21). We can simplify (35) as: 

R = (C l +C 2 + (E-E l -E 2 )&ONE)»\. (37) 
We now find the expression for E in terms of the packed 64-bit registers in (34). The 
simplification in (35), amounts to the following: 

(P+0»3 = ((P»2) + (Q»2) + £)»!, (38) 
where P and Q are unsigned integers. Let p 0 be the least significant bit of P and p x the 
next significant bit of P. Similarly, let q 0 be the least significant bit of Q and q x the next 
significant bit of Q. The simplification in (38) results in an error E when the last 2 bits of 
P and Q add up to a number > 4. The condition that determines this error E is: 

ip\ &?i)I(po&4o & (pi ki))- 

We can prove that p x , p Q , q x , q 0 can be expressed in terms of the registers in (34) as the 
least significant bits of the following packed 64-bit registers respectively: 

P X =(EC X *~(EB X \EB 2 )), 
P 0 = (EB X *EB 2 ), 

Q X =(EC 2 -~(EB 3 \EB 4 )), 
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Q 0 = (EB 3 *EB 4 ), (39) 
From (39), we can express E as the least significant bit of: 

£ = (/>!& Gi) I (A) & Qo & 1 (40) 

We note that E e {0, 1}, and£ r = (£-£, -£ 2 ) G M 0, !}• From (37), we have: 



R = 



PA VG{C X , C 2 ) - GME when £ r = -1 

PAVG(C U C 2 ) - (C, A C 2 ) & CW£ when E T = 0 . (41) 
PAVG{C X ,C 2 ) when £ r =1 



We simplify (41) as: 

* = P4 KG(Cj ,C 2 ) - ((C, A C 2 ) | (£, A £ 2 A £)) & (£, | £ 2 | ~£) & CW£, (42) 
where £, = £C, & (jE»i | EB 2 ) and £ 2 = EC 2 & (££ 3 | £5 4 ). We can further simplify (42) 
as an expression in terms of EC], EC 2 , EB\, EB 2 , £5 3 , and EB 4 so that we can skip the 
computations of E x , E 2 , and E as follows: 

U=EC l \EC 2 , 
V = EB ] \EB 2 , 
W =EB i \EB 4 , 
X = V\W, 
Y = U\X, 
Z = (EC ] &EC 2 &X), 
T = U &V&W& ((EB X & EB 2 ) | (EB 3 & EB A )), 
R=PA VG{C X ,C 2 ) - ((C, A C 2 ) \Z\T)&Y& ONE, (43) 



[58] The solution in (43) is shown in pseudo-code in Table III, below. Any suitable 
language, coding technique, circuitry or combination of hardware and software can be 
used to achieve the functionality shown in the pseudo-code presented herein. The 
approach of Table III uses 32 instructions as compared to the conventional 35 
instructions. The count of 32 instructions is obtained by counting each logical and 
arithmetic operation of (43) along with those of (34). Other instruction counts in this 
application are obtained, similarly. Clearly, the approach of Table III is not as efficient 
as the Type 2 algorithms. However, there are at least 2 benefits of this approach: 
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(1) We can systematically arrive at approximate solutions by making assumptions 
on the error/correction terms EC h EC 2 , EB h EB 2 , EB Z , and EB 4 (see Section 
3.A.iii). 

(2) In special cases, where various Afs are same, we can simplify the computation 
considerably and obtain efficient exact and approximate solutions (see Sections 
3.F-3.H). 



#define P(a,b) (((a) + (b) + 1) » 1) 

M=P(a01,a02); 

b2 = P(a03,a04); 

b3 = P(a05,a06); 

b4 = P(a07,a08); 

cl=P(bl ,b2 ); 

c2 = P(b3 ,b4 ); 

d =P(cl ,c2); 

ebl = aOl A a02; 

eb2 = a03 A a04; 

eb3 = a05 A a06; 

eb4 = a07 A a08; 

ecl=bl A b2; 

ec2 = b3 A b4; 

ed =cl A c2; 

u = eel | ec2; 

v = ebl | eb2; 

w = eb3 | eb4; 

x = v | w; 

y =u|x; 

z = eel & ec2 &x; 

t = u & v & w & ((ebl & eb2) | (eb3 & eb4)); 
e = ((ed & y) | z 1 1) & 0x01 ; // Exact solution 
xl =CLIP(d-e); 

TABLE III 



iii. Approximate SIMP Solution 

[59] We have many approximate solutions by assuming the least significant bit of EB\ , 
EB 2 , EBi, EB4, EQ, or EC2 as 0 or 1. With the assumption that the least significant bit of 
EBi = 1, and EB 2 = EB 3 = EB 4 = 0, we get from (43): 

R = PA VG(C X ,C 2 ) - ((C, A C 2 ) I (£C, & EC 2 )) & ONE. (44) 
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[60] An example pseudo-code implementation of this solution is shown in Table IV, 
below. This approach uses 14 instructions, and produces a maximum error of ±1 in the 
final result for 9.38% of all possible values of A x ,...^ig between [0,255]. The error never 
exceeds ±1. This solution with 14 instructions, and a maximum error of ±1 for less than 
l/10 th of the data is acceptable in many applications like post-processing, where a 
difference of 1 gray value in the displayed frame is imperceptible to most of us. Yet, we 
receive a computational advantage of 35:14. 

[61] The second approximate solution makes the assumption EB\ = 0, and EBi = 1 . It 
produces the following solution: 

T = (EC X | EC 2 ) & (EB 3 | EB 4 ) & EB 2 & EB 4 , 
R = PA VG(C { ,C 2 ) - ((C, A C 2 ) | {EC X & EC 2 ) \ T) & ONE, (45) 
This solution requires 22 instructions, and produces a maximum error of ±1 in the final 
result for 6.25% of all possible values oiA h ...As between [0,255]. 



#define P(a,b) (((a) + (b) + 1 ) » 1 ) 

bl=P(a01,a02); 

b2 = P(a03,a04); 

b3=P(a05,a06); 

b4 = P(a07,a08); 

cl =P(bl ,b2 ); 

c2 = P(b3 ,b4 ); 

d=P(cl,c2); 

ecl=bl A b2; 

ec2 = b3 A b4; 

ed = cl A c2; 

e = (ed | (eel & ec2)) & 0x01 ; //approx - 9.375% 
xl =CLIP(d-e); 

TABLE IV 



B. Type 3. Filter 2: R= (A i + A , + A,+ A, + Ae + A^ + A »r + Ao + 3*ONE) » 3 
[62] We require 35 instructions to compute this filter by conventional SIMD methods. 
For the new SIMD solution, we write the filter as: 

[63] R = (((^ 1 +y4 2 4v4 3 +^4+2*OA^»2) + ((A 5 +A 6 +A 7 +A s +ONE)»2) + E)»l, 
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where E is the error/correction term that is necessary for dividing up the expression into 
two parts each containing a Type 2 Filter. We have: 

R = (C, + C 2 + (E - E x - E 2 )&ONE) » 1, 
(^ 1 +^ 2 +^3+^ 4 +2*OAr£)»2 = C x - (£, & ONE), 

{A 5 +A 6 +A 7 +A s +ONE)»2 = C 2 - (E 2 & ONE), 
E x = EC X & {EB X | EB 2 ), E 2 = £C 2 | {EB 3 & £B 4 ), 
P, = {EC X A ~{EB X | EB 2 )), P 0 = (EB X A EB 2 ), 
Q\ = (^Q ^ (EB 3 & £5 4 )), 0 O = ~(EB 3 A £5 4 ), 

£=(p 1 &e I )i(^o&eo&(^iiei))- 

Here E x and £ 2 are error/correction terms obtained from (21) and (27) respectively. 
Defining E T =(E-E X - E 2 ) e {-2, -1,0}, we have: 

PA VG(C X , C 2 ) - (C, A C 2 ) & 0A>E when E T = 0 
= ] PAVG(C x ,C 2 )-ONE when£ r =-l. (46) 

iM FG(C, , C 2 ) - CW£ - (C, A C 2 ) & CW£ when E T = -2 

We simplify (46) as: 

S=(E X *E 2 *E), 

R = PA VG{C X ,C 2 ) - {{E x &E 2 &~E)\S)& ONE - ((C, A C 2 ) & ~S) & ONE, (47) 
This solution can be further simplified as: 

P = EB 3 &EB 4 , 

U=EB X &EB 2 &P, 
V=EC X &EC 2 , 
W=EB 3 \EB A , 
X= {EC X | EC 2 ) & {{EB X & {EB 2 \ W)) \ (EB 2 &W)\P), 
Y={X\V\U), 
ED = {C X *C 2 ), 

R = PA VG{C X ,C 2 ) - {{ED | Y) & ONE) -{U&V&ED& ONE). (48) 
The solution in (47) requires 35 instructions, same as the conventional 35 instructions. 
[64] An approximate solution of (47) can be obtained with the assumption that the 
least significant bit of EB X = EB 2 = 1, and EB 3 = EB 4 = 0 is: 
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R = PA VG\C X ,C 2 ) - ((C, A C 2 ) | £C, | £C 2 ) & OM (49) 
[65] This solution requires 14 instructions, and produces a maximum error of ±1 in the 
final result for 9.38% of all possible values ofA x ,...^L$ between [0,255]. We receive a 
computational advantage of 35:14. 

[66] The second approximate solution makes the assumption EB1 = 1 , and EB2 = 0. It 
produces the following solution: 

Y = ((£C, | EC 2 ) & (EB 3 | EB 4 )) \ (EC X & EC 2 ), 
R=PA VG(C X ,C 2 ) - ((C, A C 2 )\Y)& 0W£. (50) 
This solution requires 20 instructions, and produces a maximum error of ±1 in the final 
result for 6.25% of all possible values of A h ...^4% between [0,255]. 

C. Type 3, Filter 3: R = (A , + A 2 + + A, + A c + A* ± A , + A« + 2*ONE) » 3 
[67] We require 35 instructions to compute this filter by conventional SBV1D methods. 
For the new SIMD solution, we write the filter as: 

R = ((iA x +A 2 +A 3 +A 4 +ONE)»2) + ((A 5 +A 6 +A 7 +A 8 +ONE)»2) +E)» 1, 
where E is the error/correction term that is necessary for dividing up the expression into 
two parts each containing a Type 2, Filter 2 (27). We have: 
[68] R = (C, + C 2 + (E - E x - E 2 )&ONE) » 1 , 

(A x +A 2 +A-!+A A +ONE)»2 = C X - (E x & ONE), 
(A 5 +A 6 +A 7 +A s +ONE)»2 = C 2 - (E 2 & ONE), 
E x = EC X | (EB X & EB 2 ), E 2 = EC 2 \ {EB l & EB A ), 
P, = (£C, A (EB X & EB 2 )), P 0 = ~(EB X A EB 2 ), 
Q x = {EC 2 A (EB 3 & EB 4 )), Q 0 = ~(EB 3 A ^ 4 ), 
^ = (/>, & G,) | (/> 0 & Q 0 & (^i I Qi))- 

Here E x and jfs 2 are error/correction terms obtained from (27). Defining E T =(E-E X - 
E 2 )e {-2, -1, 0}, we have the same expression for R as in (46), which we simplify as: 

S=(E X *E 2 ), 

R = PAVG\C X ,C 2 ) - HE &S)&(E X \E 2 \E))& ONE - ((C, A C 2 ) & ~(5 A E)) & ONE. 

(51) 
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[69] This solution can be further simplified as: 

P = EB X \EB 4) 
Q = EB 3 \EB 2 , 
U = (EB 2 & EB 3 & P) | (EB 4 & EB X & Q), 
V=EC X &EC 2 , 

W=P\Q, 
ED = {C^C 2 ), 

R = PAVG(C h C 2 )-(ED\U\ V\({EC ] \ EC 2 ) & W)) & ONE - ED & U & V&ONE. 

(52) 

The solution in (52) requires 34 instructions, close to the conventional 35 instructions. 
[70] The approximate solution requires the assumption that the least significant bit of 
EB] = 1, and EB 2 = EB3 =EB 4 = 0 is: 

R = PA VG{C { ,C 2 ) - ((C, A C 2 ) I EC X \ EC 2 ) & ONE. (53) 

This solution requires 14 instructions, and produces a maximum error of ±1 in the final 
result for 9.38% of all possible values of A { ,...A% between [0,255]. We receive a 
computational advantage of 35:14. 

[71] The second approximate solution makes the assumption EB\ = 1 , and EB 2 = 0. It 
produces the following solution: 

U=EB 2 &EB 4 , 
ED = (C, - C 2 ), 

R = PA KG(C, ,C 2 ) -(ED\U\ EC X \ EC 2 ) & ONE -ED&U&EC X & EC 2 & ONE. 

(54) 

This solution requires 23 instructions, and produces a maximum error of +1 in the final 
result for 6.25% of all possible values of A x ,...^i s between [0,255]. 

D. Type 3. Filter 4; R = (A^_ + A-> + A, + A. + A e + A* + A , + Aq + ONE) » 3 

[72] We require 35 instructions to compute this filter by conventional SIMD methods. 

For the new SIMD solution, we write the filter as: 

R = (((A l +A 2 +A 3 +A 4 +ONE)»2) + ((^ 5 +J 6 +^ 7 +^ 8 )»2) +E)»\, 
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where E is the error/correction term that is necessary for dividing up the expression into 
two parts each containing a Type 2, Filter. We have: 

R = (C, + C 2 + {E-E x -E 2 )&ONE) » 1, 
(A x +A 2 +A 3 +A 4 +ONE)»2 = C, - (E x & ONE), 

(A 5 +A 6 +A 7 +A s )»2 = C 2 - (E 2 & ONE), 
E x = £Ci | {EB X & £5 2 ), £ 2 = EC 2 \ EB l \ EB 4 , 
P x = (EC X A (£5, & EB 2 )), P Q = ~{EB X A EB 2 ), 
Q x = (£C 2 A (£5 3 | £5 4 )), g 0 = (EB 3 A £fi 4 ), 

£:= (Pi & qo i (/> 0 & Q 0 & (^i I <2i)). 

Here E\ and Ei are error/correction terms obtained from (27) and (30) respectively. 
Defining ET= (E-E\- E 2 )e {-2, -1,0}, we have the same expression for R as in (46), 
which we simplify as: 

S=(E 2 *E), 

R=PA VG(C X ,C 2 ) -(E X \S)& ONE - ((C, A C 2 ) & ~(£, A S)) & ONE, (55) 
This solution can be further simplified as: 

P = EB 3 &EB 4 , 
Q = EB 2 \EB 4 , 
U=(EB X & (EB 2 | 0) | (EB 2 &Q)\P, 
V=EB X &EB 2 &P, 
W=EC X \EC 2 , 
Z=(EC X &EC 2 &U), 
ED = (C X *C 2 ), 

R=PAVG(C X ,C 2 )-(ED\U\ W) & ONE - ED & {{W & V)\ Z) & ONE. (56) 
The solution in (56) requires 35 instructions, same as the conventional 35 instructions. 
The approximate solution requires the assumption that the least significant bit of EB\ = 
EBi = 0, and £5 3 = EB 4 =l is: 

R = PA VG(C X ,C 2 ) - ONE - ((C, A C 2 ) & EC X & EC 2 ) & ONE. (57) 
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This solution requires 15 instructions, and produces a maximum error of ±1 in the final 
result for 9.38% of all possible values of A x ,...^i% between [0,255]. We receive a 
computational advantage of 35:15. 

[73] The second approximate solution makes the assumption EB\ = 1, and EB 2 = 0. It 
produces the following solution: 

Q = EB 3 \EB A , 
Z = (EC X &EC 2 &Q), 
ED = (C X *C 2 ), 

R = PA VG(C X ,C 2 ) - ((ED | g | £C, | EC 2 ) & ONE) -(ED&Z& ONE). (58) 
This solution requires 23 instructions, and produces a maximum error of ±1 in the final 
result for 6.25% of all possible values ofA x ,...^ between [0,255]. 

E. Type 3, Filter 5; R = (A , + A, + A-, + A 1 + Ac + A* + A t + A«) » 3 

[74] We require 33 instructions to compute this filter by conventional SDVID methods. 
For the new SIMD solution, we write the filter as: 

R = (((A x +A 2 +A 3 +A 4 )»2) + ((A 5 +A 6 +A 7 +A 8 )»2) + E)»1, 

where E is the error/correction term that is necessary for dividing up the expression into 
two parts each containing a Type 2, Filter (30). We have: 

R = (C X + C 2 + (E-E X - E 2 )&ONE) » 1, 

(A x +A 2 +A 3 +A 4 )»2 = C, - (E x & ONE), 
(A 5 +A 6 +A 7 +A 8 )»2 = C 2 - (E 2 & ONE), 
E x = EC X | EB X | EB 2 , E 2 = EC 2 \ EB 3 \ EB 4 , 
P x = (EC X a (EB X | EB 2 )), P 0 = (EB X A EB 2 ), 
Q x = (EC 2 a (EB 3 1 EB 4 )), Q 0 = (EB 3 - EB 4 ), 
E = (P X &Q X )\(P 0 &Q Q &(P X \Q X )). 
Here E\ and E 2 are error/correction terms obtained from (27) and (30) respectively. 
Defining ET=(E-E X - E 2 )e {-2, -1,0}, we have the same expression for R as in (46), 
which we simplify as: 

R = PA VG(C X ,C 2 ) - (E x | E 2 ) & ONE - ((C x A C 2 ) & ~(E X *E 2 A E))& ONE, (59) 
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This solution can be simplified as: 

P = EB X \EB 4 , 
Q = EB 3 | EB 2 , 
U=P\Q, 

V = {EB 2 & EB 3 & P) | {EB A & EB X & Q), 
W=EC ] \EC 2 , 
Z = {EC X &EC 2 &U), 
ED = {C^C 2 ), 

R = PAVG{C U C 2 ) - {{ED \U\W)& ONE) - {ED & {{W& V)\Z)& ONE). (60) 
The solution in (59) requires 34 instructions, close to the conventional 35 instructions. 
[75) The approximate solution requires the assumption that the least significant bit of 
EB\ = 1, and EB 2 = EB 3 = EB 4 = 0 is: 

R = PA KG(C, ,C 2 ) - ONE - (C, A C 2 ) & EC X & EC 2 & ONE. (6 1 ) 

This solution requires 15 instructions, and produces a maximum error of ±1 in the final 
result for 9.38% of all possible values of^ lv ..^g between [0,255]. We receive a 
computational advantage of 33:15. 

[76] The second approximate solution makes the assumption EB1 = 1, and EB2 = 0. It 
produces the following solution: 

W=EC l \EC 2 , 

R = PA VG{C X ,C 2 ) - ONE - {{C x A C 2 ) & (( W & EB 3 & EB 4 ) | {EC X & EC 2 )) & ONE). 

(62) 

This solution requires 21 instructions, and produces a maximum error of ±1 in the final 
result for 6.25% of all possible values of A X ,...A% between [0,255]. 

F. Type 3. Special Filter 1 : R = ( A , + 2A, + 2A g + 2A 7 + A ? + 4*ONE) » 3 

[77] This filter is an important loop filter for de-blocking in JVT video compression 

standards. 
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i. Conventional SIMP Solution 

[78] This filter can be implemented in SMD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (5 Instructions) A iL = Unpack Low 4 Bytes of A h for i e { 1 ,2,3,5,7} , 

2. (5 Instructions) A iH = Unpack High 4 Bytes of A t , for i e {1,2,3,5,7}, 

3. (9 Instructions) Add and Shift lower 4 words of A],...,A S to obtain lower 4 words of 
i? L as: 

R L = (A iL + 2A 3L + 2A 5L + 2A 1L + A 2L + A*ONE A ) » 3, 

4. (9 Instructions) Add and Shift higher 4 words of A } ,...J. 5 to obtain higher 4 words of 
i? H as: 

R u = (A lH + 2A m + 2A 5ii + 2A m + A m + 4*ONE 4 ) » 3, 

5. (1 Instruction) Pack R u and R L into final register R. 

We require 29 instructions to compute this filter by conventional SIMD methods. 

ii. Efficient SIMP Solution 
[79] From (34), we get the: 

5, = PA VG(A j A 2 ), B 2 = A 3 , B 3 =A 5 ,B 4 = A 7 , 

C, =PAVG(.B i A 3 ), C 2 = PAVG(A 5 4 7 ), 
EB ] = {AfA 2 ), EB 2 = 0, EB 3 = 0, EB 4 = 0, 
^C 1 =(5,^ 3 ), J BC 2 = (^ 7 ). (63) 

In (43) we get: 

R = PAVG(C h C 2 ) - ((Cj A C 2 ) | (£C, & £C 2 & £5,)) & (EC l \ EC 2 \ EB X ) & ONE. 

(64) 

The solution in (64) requires 16 instructions with a computational benefit of 29:16. 

iii. Approximate SIMP Solution 

[80] The approximate solution with the assumptions that the least significant bit of EB\ 
= EC\ = 0, and ECi = I is: 

R = PA FG(C„C 2 ) - (C, A C 2 ) & ONE. (65) 
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It requires 7 instructions, and produces a maximum error of ±1 in the final result for 
12.5% of all possible values of A U ...,A S between [0,255]. The computational advantage is 
29:7 (approx. 4 times speedup). 

G. Type 3, Special Filter 2: R = (A i + A-> + A? + 3A,« + 2A i + 4*ONE) » 3 

[81] This filter is also an important loop filter for de-blocking in the JVT video 
compression standard. 

i. Conventional SIMP Solution 

[82] This filter can be implemented in SIMD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (5 Instructions) A iL = Unpack Low 4 Bytes of A i9 for i e { 1 ,2,3,4,7} , 

2. (5 Instructions) A iH = Unpack High 4 Bytes of A i9 for i e { 1 ,2,3,4,7} , 

3. (8 Instructions) Add and Shift lower 4 words of ^ 1? ...^4 5 to obtain lower 4 words of 
i? L as: 

*L = iA\L + ^2L + ^3L + 3^4L + ^7L + 4*ONE 4 ) » 3, 

4. (8 Instructions) Add and Shift higher 4 words of ^ lv ..^4 5 to obtain higher 4 words of 
R u as: 

R H = (A 1H + v4 2H + ^ 3 H + 3^4H + 2^ 7 H + 4* ONE J » 3, 

5. (1 Instruction) Pack R H and 7? L into final register R. 

[83] We require 27 instructions (including two multiplications by 3) to compute this 
filter by conventional SIMD methods. 

ii. Efficient SIMD Solution 
[84] From (34), we get the: 

B x = PAVG(A l9 A 2 ) 9 B 2 = PAVG(A 39 A 4 ), B 3 = A 4 , B 4 = A 7 , 
Cj = PAVG(B l9 B 2 ), C 2 = PAVG(A 49 A 7 ), 
EB X = (A { A A 2 ), EB 2 = (A 3 A A 4 ), EB 3 = 0, EB 4 = 0, 

£C, = (5^), £C 2 = (^ 4 ^ 7 ). (66) 

In (43) we get: 
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S=(EB l \EB 2 ) 9 

R = PA VG{C X ,C 2 ) - ((Q A C 2 ) | (£Q & EC 2 &S))& (EC X \EC 2 \S)& ONE. (67) 
The solution in (65) requires 19 instructions with a computational benefit of 27:19. 

iii. Approximate SIMP Solution 

[85] The first approximate solution with the assumptions that the least significant bit of 
EC X =EC 2 = 1, and EB X = EB 2 = 0 is: 

R = PA VG(C h C 2 ) - (Cj A C 2 ) & ONE. (68) 
It requires 8 instructions, and produces a maximum error of ±1 in the final result for 
12.5% of all possible values of ^ lv ..y4 8 between [0,255]. The computational advantage is 
27:8. 

The next approximate solution is with the assumption that the last bit of EB { = EB 2 = 
1, which gives us: 

R = PAVG{C h C 2 )-~({CfC 2 ) | (EC { &EC 2 )) & ONE. (69) 
This solution requires 12 instructions and produces a maximum error of ±1 in the final 
result for 6.25% of all possible values of A\ 9 ... 9 A% between [0,255]. The computational 
advantage is 27:12. 

H. Type 3, Special Filter 3: R = (A , + A 2 + A, + 2A A + A g + A^ + A ^ + 4*ONE) » 3 
[86] This filter is used for de-blocking in post-processing. 

i. Conventional SIMP Solution 

[87] This filter can be implemented in SMD architecture (assuming sufficient 
memory) by using the following steps: 

I. (7 Instructions) A iL = Unpack Low 4 Bytes ofA i9 for i € {1,2,3,4,5,67}, 

2. (7 Instructions) A iH = Unpack High 4 Bytes of A i9 for i e {1,2,3,4,5,67}, 

3. (9 Instructions) Add and Shift lower 4 words of ^ l5 ...^ 5 to obtain lower 4 words of 
R L as: 

R L = (4l L + *2L + ^3L + 2^4L + ^5L + ^6L + ^7L + 4*ONE 4 ) » 3, 
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4. (9 Instructions) Add and Shift higher 4 words of A } ,...y4 5 to obtain higher 4 words of 
R H as: 

R H = (A ]H + A m + A m + 2A 4H + A 5H + A 6H + A 1H + A*ONE 4 ) » 3, 

5. (1 Instruction) Pack R u and R L into final register R. 

We require 33 instructions to compute this filter by conventional SIMD methods. 

ii. Efficient SIMP Solution 
[88] From (34), we get the: 

fl, = PA VG(A , J 2 ), B 2 = PA VG(A 3 J 5 ), B 3 = PA VG(A 6 J 7 ), B A = A 4 , 

C, = PAVG(B U B 2 ), C 2 = PAVGiB i 4 4 ), 
EB { = (A^A 2 ), EB 2 = (A 3 *A 4 ), EB 3 = {A^A n ), EB 4 = 0, 

EC X = (B^B 2 ), EC 2 = (BfA 4 ). (70) 

In (43) we get: 

U=EC X \EC 2 , 
V=EB X \EB 2 , 
X = V\EB lt 
Y =U\X, 
Z =(EC { &EC 2 &X), 
T =U &V&EB l &EB 2 &EB i , 
R = PA VG{C X ,C 2 ) - ((C, A C 2 ) | Z | T) & Y & ONE, (7 1 ) 

The solution in (71) requires 27 instructions with a computational benefit of 33:27. 

iii. Approximate SIMD Solution 

[89] We get an approximate solution with the assumptions that the least significant bit 
of EB 2 = 1, and EB 2 = EB 3 = 0 as: 

R = PA VG(C X ,C 2 ) - ((C, A C 2 ) | (£C, & EC 2 )) & ONE. (72) 

The solution in (72) requires 13 instructions, and produces a maximum error of ±1 in the 
final result for 6.25% of all possible values of A h ...yA-i between [0,255]. The 
computational advantage is 33:13. 
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I. Type 3, Summary of Results 

[90] Table V below summarizes the instructions required to compute each filter (given 
sufficient memory) by the efficient and conventional SMD methods. For the Efficient 
case, we give the instructions required for the exact and approximate solutions. 



Summary of Results for Type 3 FIR Filters. 



Type 3 Filters 


Conventional 
Method 


Efficient Method 


Speedup 


Exact 


Approx. 


Exact 


Approx. 


{A x + A 2 + ... + A s +4*ONE)»3 


35 


32 


14 


1.1 


2.5 ; 


(A x +A 2 + ...+A 8 + 3*ONE)»3 


35 


35 


i4 


1.0 


2.5 :■ 


(A { +A 2 + ... +Ag + 2*ONE)»3 


35 


34 


: 14 


1.0 


• . 2.5 . ; 


(A { +A 2 + ... +A S + \*ONE) » 3 


35 


35 


'i ■■" ; i-5;. '■■< 


1.0 


f -2.3,' :,; 


(A l +A 2 + ...+A s )»3 


33 


34 


U:-: ; ri5."|;J 


1.0 


2.2 


(A l +2A 3 +2A 5 +2A 7 +A 2 +4*ONE)»3 


29 


; 16 




1.8 


. 4.1 . 


(A l +A 2 +A 3 +3A 4 +2A 7 +4*ONE)»3 


27 


.19 • 




1.4 


3.4 


(A [ +A 2 +A 3 +2A 4 +A 5 +A 6 +A 7 +4*ONE)»3 


33 


27 


13 , 




2.5 :: 



TABLE V 

[91] The shaded areas show significant improvements in efficiency due to the analyses 
developed here. 



4. Type 4 FIR Filters 

[92] There are 9 different Type 4 FIR filters depending on the 9 choices of c in (9). For 
the sake of brevity, we shall only discuss the case of c=8. We define the following packed 
64-bit registers, each containing 8 data elements of one byte each: 

B l =PAVG(A^ 2 ),B 2 = PAVG(A^ A ), B 3 = PAVG(A 5 ^ 6 ), B 4 = PAVG(A 7 J S ), 

B 5 = PA VG{A 9t A xo ), B 6 = PA VG{A n A X2 ), B 7 = PAVG{A {3 4\^ B s = PA VG(A l5 4i 6 ), 
C, = PA VG(B { ^ 2 ), C 2 = PA VG{B^), C 3 = PA VG(B 5 ,B 6 ), C 4 = PA VG{B 7 £z), 
D x =PAVG(C h C 2 ),D 2 = PAVG(C 3 ,C 4 ), 
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EB X = (AfA 2 ), EB 2 = (AfA 4 \ EB 3 = (A 5 A A 6 ), EB 4 = (AfA s ), 
EB 5 = (A 9 *A X0 ), EB 6 = (A x X *A X2 ), EB 1 = (A X ^A X4 ), EB Z = (A X5 *A X6 ), 
EC X = (B X *B 2 ), EC 2 = (BfB 4 ), EC 3 = (B 5 *B 6 ), EC 4 = (P 7 *5 8 ), 
ED X =(C X ^C 2 \ED 2 = (C^C 4 ), 

E x = EC i & (EB X | EB 2 ), E 2 = EC 2 & (EB 3 | EB 4 ), 
E 3 = PC 3 & (EB 5 | EB 6 ), E 4 = EC 4 & (EB 7 \ EB S ), 

P, = (EC X * ~(£P, | £5 2 )), P 0 = (EB X * EB 2 ), 
Q x = (EC 2 * ~(PP 3 1 P5 4 )), Q 0 = (EB 3 - EB 4 ), 
ER X = (P, & Q x ) | ((P, | 0,) & P 0 & Q 0 ), 

R x = (EC 3 * ~(EB 5 1 P5 6 )), i? 0 - {EB 5 * EB 6 ), 
S x = (EC 4 A ~(£5 7 | ^g)), S 0 = (£5 7 - £5 8 ), 
ER 2 = (R X &S X )\((R X \S X )&R 0 &S 0 ), 

C/ 2 = £2), * £, * £ 2 a ((P, & e ,) | ((P, |2i)&P 0 & e 0 )), 
C/^P^e.-CPo&go), 

V 2 = ED 2 -E^E 4 - ((R x & S x ) | ((/?, | S x ) & i? 0 & S 0 ))> 

E = (U 2 & V 2 ) | (((7 2 | F 2 ) & £/, & F,) | ((t/ 2 | r 2 ) & ((7, | V x ) &U 0 & Fq), 

£P, = (ED X | (£, A P 2 A £»,)) & C^i I El I 
ZTP 2 = {ED 2 | (£ 3 ^ £ 4 * £Rj)) & (£ 3 | E 4 | ~£tf 2 ). 
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A. Type 4, Filter 1: R = ( A } ± A , + ... + A i g + A - ^ + 8*ONE) » 4 
1 Conventional SIMP Solution 

[93] This filter can be implemented in SMD architecture (assuming sufficient 
memory) by using the following steps: 

1. (16 Instructions) A iL = Unpack Low 4 Bytes of A i9 for i = 1,. . .,16, 

2. (16 Instructions) A iH = Unpack High 4 Bytes of A i9 for i = 1,. . .,16, 

3. (17 Instructions) Add and Shift lower 4 words of A x , . . . >A x 6 to obtain lower 4 words of 
R L as: 

tf L = (A XL + ^ 2L +...+^ 15L + ^ 16L + 4*ONE 4 ) » 4, 

4. (17 Instructions) Add and Shift higher 4 words of ^ lv ..^4 16 to obtain higher 4 words 
of Z? H as: 

R H = (A lH + ^ 2H + . . . + A l5H + ,4 16H + 4*0^ 4 ) » 4, 

5. (1 Instruction) Pack R H and i? L into final register R. 

We require 67 instructions to compute this filter by conventional SMD methods. 
ii. Efficient SIMP Solution 

[94] In order to implement this filter efficiently, we simplify it as follows: 
R = {{A x +A 2 +. . .+^ 7 +^ 8 +4*OAT£)»3 + (A 9 +A x0 +. . .+A ls +A l6 +4*ONE)»3 + E) » 1 , 

(74) 

where £ is the error/correction term that is necessary for dividing up the expression into 
two parts each containing a Type 3, Filter 1 (42). Note that according to (42) and the 
registers in (73), we have: 

{A x +A 2 +. . .+v4 7 +4 8 +4*CW£)>>3 = D x - (£7^ & CW£), 
(i4 9 +i4 10 +...+ i 4 15 + i 4 16 +4*OJVjB)»3 = Z) 2 ~ ( ET 2 & 0A/E )> ( 75 ) 
where ET X and £7 2 are error/correction terms obtained in (42). We can simplify (74) as: 

R = (D x + Z> 2 + (E- ET X - ET 2 )8cONE) » 1 . (76) 
Note that the expressions for E, ET X , and ET 2 are given in (73). We note that E, ET X , ET 2 
€ {0, 1 } , and Ej = (is — ET X — ET 2 ) e {-2,-1, 0, 1}. From (76), we have: 
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R = 



(77) 



(78) 



PAVG(D h D 2 )-ONE-(D ] A D 2 )&ONE when E T =-2 
PAVG(D u D 2 )-ONE when £ r =-1 

PA VG(D l ,D 2 ) - (Z>, A D 2 ) & CW£ when £ r = 0 ' 
^(Z), , Z> 2 ) when £ r = 1 

We simplify (77) as: 

R = PA VG{D X J) 2 ) - {{ET X & ET 2 ) \ ~E) & {ET X \ET 2 \E)& ONE 
- {D X A D 2 ) & ~{ET { A ET 2 A E) & ONE. 
We can further simplify (78) as: 

£/, = At least 1 ED, 
U 2 = At least 1 EC, 
t/ 3 = At least 1 EB, 

{y 4 = Both£Z)s, 
U s = At least 2 ECs, 
U 6 = At least 3 ECs, 

U 7 = Ml4ECs, 
U 8 = At least 3 EBs, 
U 9 = At least 5 EBs, 
£/ 10 = Atleast7 EBs, 
U n =At least 1 ED, EC or EB, 
E\ = (EX& U u ) | (£/ 4 & (C/ 2 | t/ 3 )) | (C/, & l/ 6 ) | (C/, & £/ 5 & C/ 3 ) | (C7, & t/ 2 & U 8 ) \ 
(U,&U 9 )\(U 7 &U 3 )\(U 6 &U B )\(U 5 &U 9 )\(U 2 &U l0 % 
E 2 = (EX& U 4 & ((U 7 & U 3 ) | (U 6 & U 8 ) | (U 5 & £/ 9 ) | (U 2 & U l0 ))) | 
t/j & ((t/ 7 & t/ 9 ) | (C/ 6 & C/ 10 )), 

E=E X +E 2 . (79) 
[95] Clearly, (79) is an inefficient solution and is useful in special cases and for 
approximate solutions only. 
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iii. Approximate SIMP Solution 

[96] We have many approximate solutions by assuming the least significant bit of EB U 
EB Z , EC U ... EC 4 , ED h or ED 2 as 0 or 1. With the assumption the last bit of E = ET X = 
ET 2 = l,we get from (78): 

R = PA VG(D X J) 2 ) - ONE. (80) 
This solution requires 16 instructions, and produces a maximum error of ±1 in the final 
result for 8.6% of all possible values of A b ...^l l6 between [0,255]. The error never 

exceeds ±1. We receive a computational advantage of 67:16, and approximate 4 times 
speedup. 

B. Type 4, Special Filter 1 : R = (A t + 4 A -> + 6A^ + 4A £ + A g + 8*ONE) » 4 
[97] This filter is a Gaussian approximation filter used for post-processing. 

i. Conventional SIMP Solution 

[98] This filter can be implemented in SIMD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (5 Instructions) A iL = Unpack Low 4 Bytes of A i9 for i e { 1 ,2,3,4,5 } , 

2. (5 Instructions) A iH = Unpack High 4 Bytes of A h for / e {1,2,3,4,5}, 

3. (9 Instructions) Add and Shift lower 4 words of A^ 9 A 5 to obtain lower 4 words of 
R L as: 

Rl = (A lL + 4A 2L + 6A 3L + 4A 4L + A 5L + %*ONE A ) » 4, 

4. (9 Instructions) Add and Shift higher 4 words of A h A 5 to obtain higher 4 words 
ofi? H as: 

R H = (A m + 4A 2H + 6A m + 4^ 4H +^ 5H + 8*OJV£ 4 ) » 4, 

5. (1 Instruction) Pack R H and R L into final register R. 

We require 29 instructions to compute this filter by conventional SIMD methods. 

ii. Efficient SIMD Solution 

From (73), we get the: 
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B X =PAVG(A X J 2 ),B 2 =A 3 , B 3 =A 3 ,B 4 = A 3 ,B 5 =A 4 ,B 6 =A 4 , B n =A s ,B % =A 5 , 
C, = PA VG(B X J 3 ), C 2 = A 3 ,C 3 = A 4 , Q = A 5 , 

D x =PAVGiC x A 3 ),D 2 = PAVG{A 4y A s \ 
EB X = A X A A 2 , EB 2 = EB 3 = EB 4 = EB 5 = EB 6 = EB 7 = EB S = 0, 
EC X = (,B X A A 3 ), EC 2 = EC 3 = EC 4 = 0, 
ED X = C X A A 3 , ED 2 = A 4 A A 5 , 

E x — EC] & EB\, E 2 — E 3 — E 4 — 0, 
P, = EC] A ~EB X , P 0 = EB X , Q x = 1, & = 0, = Pi, 
R { = \, R 0 = 0,S x = l,S 0 = 0,ER 2 = l, 
U 2 = ED X A E X A P X , U x =~P„ C/ 0 = P 0 , 
F 2 = ~£D 2 ,Fi=0,F 0 -0, 
E=U 2 &V 2> 

ET X = (ED X | (E x A £/?,)), PT 2 = 0. (81) 

From (78) we get: 

R = PA VG(D X J) 2 ) - (~E & ET X & ONE) - (D X A D 2 ) & ~{ET X A E) & ONE. (82) 
We can simplify (82) as follows: 

U=EC X \EB X , 

R = PAVG(D X J) 2 ) - (((D X A D 2 ) & (ED X \ ED 2 \ U)) \ (ED X & ED 2 &U))& ONE. 

(83) 

The solution in (83) requires 19 instructions with a 29:19 computational advantage. 
iii. Approximate SIMP Solution 

[99] We can assume the least significant bit of ED X , EDi, EC\, or EB\ as 0 or 1 to get 
several approximate solutions. We first make the assumption that the least significant bit 
of ED\=\, and ED 2 = EC\ = EB X = 0, to get the following solution: 

R = PA VG(D X J) 2 ) - (D X A D 2 ) & ONE. (84) 
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This solution requires 8 instructions, and produces a maximum error of ±1 for 12.5% of 
all possible values of A h .,. 9 A 5 between [0,255]. The computational advantage is 29:8 
(more than 3 times speedup). 

The second approximate solution makes the assumption that the least significant bit of 
EC i = 1, and EB x = 0, to get the solution: 

R = PAVG(D X JD 2 ) - ((D X A D 2 ) \ (ED ] & ED 2 )) & ONE. (85) 
This solution requires 12 instructions, and produces a maximum error of ±1 for 6.25% of 
all possible values of ^ lv ..^4 5 between [0,255]. The computational advantage is 29:12. 

C. Type 4, Special Filter 2: 

R = (A 1 +A^ +2A g +2A*+2A 7 +2A B H-4A 0 + A^+A^+8*ONE) » 4 
[100] This filter is also a Gaussian approximation filter used for post-processing. 

L Conventional SIMP Solution 

[101] This filter can be implemented in SIMD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (9 Instructions) A iL = Unpack Low 4 Bytes of A i9 for i e { 1 ,2, ... ,9} , 

2. (9 Instructions) A iH = Unpack High 4 Bytes of A i9 for i e { 1 ,2, ... ,9} , 

3. (15 Instructions) Add and Shift lower 4 words of A u A 5 to obtain lower 4 words of 
R L as: 

R L = (A XL + A 2L + 2A 5L + 2A 6L + 2A 7L + 2A SL + 4A 9L + A 3L + A 4L + S*ONEJ » 4, 

4. (15 Instructions) Add and Shift higher 4 words of A x , A 5 to obtain higher 4 words 
of 7? H as: 

i? H = 04 1H + A 2H + 2^ 5H + 2A 6H + 2^ 7H + 2^ 8H + 4A 9H + ^ 3H + ^ 4H + 8*CW£ 4 ) » 4, 

5. (1 Instruction) Pack i? H and R L into final register /?. 

We require 49 instructions to compute this filter by conventional SIMD methods. 

ii. Efficient SIMD Soluti n 
[102] From (73), we get the: 
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B X =PAVG(A X J 2 ),B 2 = PAVG(A 3 J 4 ), 2? 3 = A 5 , 2? 4 = A 6 , B s = A 7 , B 6 = A%, B n =A 9 ,B % = 

A 9 , 

C, = PA VG(B X ,B 2 ), C 2 = PA VG(A 5 A 6 ), C 3 = PA VG(A 7 ^), C 4 = A 9 , 
D x =PAVG(C h C 2 ),D 2 = PAVG(C 3 ^ 9 ), 
EB X = (A X A A 2 ), EB 2 = (A 2 A A 4 ), EB 3 = EB A = EB 5 = EB 6 = EB 1 = EB & = 0, 
EC X = (B X A B 2 ), EC 2 = (A 5 -A 6 l EC 2 = (A 7 *A 8 ), EC 4 = 0, 
ED X = (C X A C 2 ),ED 2 = (C 3 A A 9 ), 

E x = EC X & (EB X | EB 2 ), E 2 = E 3 =E A = 0, 
Pj = EC X A | £Z? 2 ), P 0 = (EB X A EB 2 ), Q x = ~EC 2 , Q 0 = 0, ER X =P X &Q X , 
R X =~EC 3 , Rq = 0,S x = 1,S 0 = 0,ER 2 =R x , 
U 2 = ED X A E x * (P x & Q x ), U X =P X A Q X ,U 0 = P 0 , 
V 2 = ED 2 A R x ,V x =~R x ,V 0 = 0, 
E = (U 2 &V 2 )\((U 2 \V 2 )&U X &V X ), 
ET X = (ED X | (E x A ER X )) & (£, | ~ER X ), ET 2 = (£D 2 & ~£/? 2 ). (86) 
[103] The final solution is same as (78). We can simplify this solution as follows: 

U=EB X \EB 2 , 

V=EC X \EC 3 , 
W=EC 2 \U\V, 
Z = ED X \ED 2 , 

F = Z\W, 
H = EC X &EC 3 , 
G = (EC 2 &V)\H, 

R = PA VG(D X J) 2 ) - (((D X A D 2 )&F) \ (ED X &ED 2 &W) | (Z & ((EC 2 &H)\(G&U)))) & ONE. 

(87) 

[104] The solution in (87) requires 36 instructions with a 49:36 computational 
advantage. 
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iii. Approximate SIMP Soluti n 

[105] We suggest 2 approximate solutions for this filter. For the first approximate 
solution, we assume the least significant bit of ECi=0 9 and EB\= EB 2 =l to get the 
following: 

R = PAVG(D ly D 2 ) - ((D, A Z> 2 ) I (ED } &ED 2 ) | ((ED X \ED 2 ) & EC X & EC 2 )) & ONE. (88) 
This solution requires 21 instructions, and produces a maximum error of ±1 for 6.25% of 
all possible values of A l9 ...,A 9 between [0,255]. The computational advantage is 49:21 
(more than 2 times speedup). 

[106] The second approximate solution makes the assumption that the least significant 
bit of EB X = EB 2 - 1. We get the solution: 

[/= (EC { & (EC 2 | EC 3 )) | (EC 2 & EC 3 ) 9 
R = PA VG(D ]y D 2 ) - ((D } A D 2 ) | (ED ] &ED 2 ) \ ((ED ] \ED 2 ) & CO) & ONE. (89) 
This solution requires 25 instructions, and produces a maximum error of ±1 for 3.12% of 
all possible values of A l9 ...,A 9 between [0,255]. The computational advantage is 49:25 
(nearly 2 times speedup). 

P. Type 4, Special Filter 3: 

R = (A i +2A ? +2A ? +2A /< -t-2A g +2A^+2A >7 +2A g +Ao+8*ONE) » 4 
[107] This filter is also used for post-processing. 

i. Conventional SIMP Solution 

[108] This filter can be implemented in SIMD architecture (assuming sufficient 
memory) by using the following steps: 

1. (9 Instructions) A iL = Unpack Low 4 Bytes of A h forz e {1,2,. ..,9}, 

2. (9 Instructions) A iH = Unpack High 4 Bytes of A i9 for ie {1,2,. . .,9}, 

3. (17 Instructions) Add and Shift lower 4 words of A x > A 5 to obtain lower 4 words of 
R L as: 

*L = (^il + 2^2L + 2A 3L + 2A 4L + 2A 5L + 2A 6L + 2A 7L + 2^ 8L +^ 9L + 8*OA^ 4 ) » 

4, 
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4. (17 Instructions) Add and Shift higher 4 words of A x , A 5 to obtain higher 4 words 
of R H as: 

R H = (A lH + 2A 2H + 2A m + 2A 4H + 2A 5H + 2A 6H + 2A m + 2A m + A 9H + S*ONE 4 ) » 

4, 

5. (1 Instruction) Pack R H and R L into final register R. 

We require 53 instructions to compute this filter by conventional SEVID methods. 

ii. Efficient SIMP Solution 
[109] From (73), we get the: 

B X =PAVG(A X J 2 ),B 2 = A 2 , B 3 = A 3 ,B 4 =A 4 ,B 5 =A 5 ,B 6 = A 6 , B 7 =A 7 ,B i =A s , 

C x = PA VG(B lf A 2 ), C 2 = PA VG(A 3 ^ 4 ), C 3 = PA VG(A 5 J 6 ), C 4 = PA VG(A 7 A%), 
D x = PA VG(C X ,C 2 ), D 2 = PA VG(C 3 ,C 4 ), 
EB X = (A X A A 9 ), EB 2 = EB 3 = EB 4 = EB 5 = EB 6 = EB 7 = EB % = 0, 
EC X = (BfA& EC 2 = (A 3 *A 4 ), EC 3 = (A 5 *A 6 ), EC 4 = (A 7 *A 8 ), 
ED X = (C,*C 2 ), ED 2 = (C 3 *Q), 

E\ — EC\ & EB\ , E2 = E^ — E^ — 0, 
P x = EC X A ~EB X , P 0 = EB X , Q x = ~EC 2 , Q 0 = 0, ER X = P x & Q x , 
i?j — ~EC^ i?o = O5 S\ = ~ECA t , Sq — 0, ER2 = R\ & S\ 9 

U 2 =ED X -E x * (P x & Q x ), U x =P X * Q x , U 0 = P 0 , 
V 2 = ED 2 * (R x & S x ), V X =R X * S x , V 0 = 0, 
E = (U 2 &V 2 )\((U 2 \V 2 )&U X &V X ), 
ET X = (ED X I (E x * ER X )) & (E x | ~ER X ), ET 2 = (ED 2 & ~ER 2 ). (90) 
The final solution is same as (78). We can simplify this solution as follows: 

U X =EC X &EC 2 , 
U 2 = EC 3 & EC 4 , 
U 3 = ED X &ED 2 , 
V X =EC X \EC 2 , 
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V 2 = EC 3 \EC 4 , 
V 2 =ED X \ED 2 , 
V 4 =V 2 \EB h 
W=(V X &V 2 & EB X ) | (U x & V 4 ) | (U 2 & {V x | EB X )\ 
F = Vi | K4, 

G =U l &U 2 &EB ] , 
H =G&U 3 &E, 

R = PAVG{D X J) 2 ) - (((£ & (K 3 I F)) I (C/ 3 & F) | (7 3 & W) \ G) & ONE) -(H& ONE). 

(91) 

[110] The solution in (91) requires 46 instructions with a 53:46 computational 
advantage. 

iii. Approximate SIMP Solution 

[111] We suggest 2 approximate solutions for this filter. For the first approximate 
solution, we assume the least significant bit of EC\ = ECi = 1, and EC?, = EC4 = 0 to get 
the following: 

R = PA VG(D l P2) ~ (Pi AD 2) I ( ED i & ED i) I i( ED i I ED 2> &EB l ))& ONE. (92) 
This solution requires 19 instructions, and produces a maximum error of ±1 for 9.38% of 
all possible values of A h ...^l 9 between [0,255]. The computational advantage is 53:19 

(more than 3 times speedup). 

The second approximate solution makes the assumption that the least significant bit of 
EC\ = l, and EC 3 = 0. We get the solution: 

W = (EC 4 & EB X ) I (EC 2 & (EC 4 \ EB X )), 

R=PA VG{D X J) 2 ) - ((£>! A D 2 ) \(ED x 8c ED 2 ) \ ({ED X \ ED 2 ) &W))& ONE. (93) 
[112] This solution requires 25 instructions, and produces a maximum error of ±1 for 

6.25% of all possible values of A\ A 9 between [0,255]. The computational advantage is 

53:25 (more than 2 times speedup). 
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E. Type 4, Special Filter 4: R = rA 1 +2A 2 +3A 2 +4A 1 +3A f +2A f +A z +8*ONE) » 4 
[113] This filter is also used for post-processing. 



i. Conventional SIMP Solution 

[114] This filter can be implemented in SEvlD architecture (assuming sufficient 
memory) by using the following steps: 

1 . (7 Instructions) A iL = Unpack Low 4 Bytes of A t , for i e { 1 ,2, ... ,7} , 

2. (7 Instructions) A iH = Unpack High 4 Bytes of A t , for i e { 1 ,2, . . . ,7 } , 

3. (13 Instructions) Add and Shift lower 4 words of A b A 5 to obtain lower 4 words of 
i? L as: 

R L = (A , L + 2A 2L + 3A 3L + U 4L + 3A 5L + 2A 6L + A 1L + S*ONE 4 ) » 4, 

4. (13 Instructions) Add and Shift higher 4 words of A x , A 5 to obtain higher 4 words 
of R u as: 

R u = (A m + 2A m + 3A m + 4A 4H + 3A 5H + 2A m + A m + 8*ONE 4 ) » 4, 

5. (1 Instruction) Pack i? H and R L into final register R. 

We require 41 instructions to compute this filter by conventional SIMD methods. 

ii. Efficient SIMP Solution 
[115] From (73), we get the: 

B x = PA VG{A , At), B 2 = A 2 , B 3 = A 4 , B 4 = A 4 , B 5 =A 3 ,B 6 =A 5 , B 7 =A 6 ,B g = 

PAVG(A 3 A 5 ), 

C, = PA VG{B X A 2 )> C 2 = A 4 , C 3 = PA VG(A 3 A 5 ), Q = PA VG{A 6t B % ), 

D x = PA VG(C { A 4 ), D 2 = PA VG(C 3 ,C 4 ), 
EB X =A X A A 7 , EB 2 = EB 3 = EB 4 = EB 5 = EB 6 = EB 7 = 0, EB S = A 3 A A 5 , 
EC X = BfA 2 , EC 2 = 0, EC 3 = A 3 A A 5 , EC 4 = A 6 A B%, 
ED X = C } A A 4 , ED 2 = C 3 A C 4 , 

E x = EC] & EB h E 2 = E 3 = 0, E 4 = EC 4 & EB B , 
P, = EC X A ~EB h P 0 = EB l ,Q x = \,Q 0 = 0, ER X = P h 
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R x = ~£C 3 , R 0 = 0, 5j = EC 4 A ~££ 8 , S 0 = EB S , ER 2 = R X &S X , 



U 2 = ED X A P,, U x = ~P„ U Q =P 0 , 
V 2 = ED 2 *E 4 *(R X &S X ),V X =R X *S X ,V 0 = S 0 , 
E=(U 2 & V 2 ) | ((U 2 | V 2 ) &U X & V x ) | (([/ 2 | V 2 ) & (f/, | V ,) & £/ 0 "& F 0 ), 
£7/, = (ED, | (E, A P,)) & (£, | ~P,), PT 2 = (PZ) 2 I (£ 4 A £P 2 )) & (£ 4 | ~ER 2 ). (94) 
The final solution is same as (78). We can simplify this solution as follows: 

U x =EC 3 \EC 4 , 
U 2 = EB X \EB S , 
U 3 =ED X \ED 2 , 
U 4 = EC X \U U 
U 5 = U 4 \U 2 , 

u 6 = u 5 \u 3 , 

U 7 =EC 3 &EC 4 , 
U 8 = (EC X &U X )\U 7 , 
U 9 = (EC X &U 7 )\(U S &U 2 ), 
R = PAVG(D X J) 2 ) - (((D X A D 2 ) & U 6 ) | (ED X & ED 2 & U 5 ) | (U 3 & t/ 9 )) & ONE. (95) 
The solution in (95) requires 36 instructions with a 41:36 computational advantage. 

iii. Approximate SIMP Solution 

[116] We suggest 2 approximate solutions for this filter. For the first approximate 
solution, we assume the least significant bit of EC\ = EC3 = 1 , and EC4 = EB\ = 0 to get 
the following: 

R = PA VG(D X J) 2 ) - ((D X A D 2 ) I (ED, & ED 2 ) \ ((££>, \ ED 2 ) & EB$) & ONE. (96) 
This solution requires 19 instructions, and produces a maximum error of +1 for 6.25% of 
all possible values of A X ,...A 7 between [0,255]. The computational advantage is 41:19 
(more than 2 times speedup). 

[117] The second approximate solution makes the assumption that the least significant 
bit of EQ = 1, and EB X = 0. We get the solution: 
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U 9 = (EC X & EC A ) | ((£C, | EC 4 ) & EB S ), 
R = PA VG(D X JD 2 ) - ((D } A D 2 ) \ (ED } & ED 2 ) \ ((ED } \ ED 2 ) & U 9 )) & ONE. (97) 
This solution requires 25 instructions, and produces a maximum error of ±1 for 3.13% of 
all possible values of A U ...,A 7 between [0,255]. The computational advantage is 41:25. 

F. Type 4, Summary of Results 

[118] Table VI below summarizes the instructions required to compute each filter 
(given sufficient memory) by the efficient and conventional SIMD methods. For the 
Efficient case, we give the instructions required for the exact and approximate solutions. 



Summary of Results for Type 4 FIR Filters. 



Type 4 Filters 


Conventional 
Method 


Efficient Method 


Speedup 


Exact 


Approx. 


Exact 


Approx. 


(A l +A 2 + ... +A l6 + S*ONE)»4 


67 


N/A 


16 ' 


N/A 


4.2 


(A , +4A 2 +6A 3 +4A 4 +A 5 +S*ONE)»4 


29 


19 


8 


1.5 


3.6 


{A ] +A 2 +2A 5 +2A 6 +2A 1 +2A s +4A 9 +A i +A /> +»*ONEy»4 


49 


36 . 


21 


1.4 


2.3 


{.A l +2A 2 +2A 1 +2A^2A 5 +2A 6 +2A 1 +2A^A g ^*ONE)»A 


53 


46 


19 


. 1.2 


2.8 


(A l +2A 2 +3A i +4A 4 +3A s +2A ( +A 1 +%*ONE)» 4 


41 


36 


19 


1.1 


2.2 



TABLE VI 



[119] The shaded areas show significant improvements in efficiency due to the analyses 
developed here. 

5. Corrected Approximation Method 

[120] Another approach, referred to as the "corrected approximation" method, can be 
used to improve upon the approaches described in sections 3 and 4 when exact 
computation is required. Steps of the corrected approximation method include the 
following: 

1 . Use PA VG to quickly get an approximation to the correct result. 

2. Adjust the approximation so that it becomes no less than the correct result. 
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3. Perform another computation to determine the value of the least significant bits of 
the correct result. 

4. Use the value of the least significant bits from step 3 to determine the error of the 
value obtained from step 2. 

5. Subtract the error value in step 4 from the approximate value in step 2 in order to 
get the correct final result. 

[121] Examples of the corrected approximation method are given in the following 
subsections for type 3 and type 4 filters, respectively. 

A. Corrected Approximation Method Applied to Type 3 Filters 

[122] As in section 3, we define the following packed 64-bit registers, each containing 8 

data elements of one byte each: 

B, =PAVG(A l9 A 2 lB 2 = PAVG(A 39 A 4 ), B 3 = PAVG(A 59 A 6 ), B 4 = PAVG(A 7 ^), 

C, = PA VG(B { ,B 2 \ C 2 = PA VG(B 3 ,B 4 ), D = PA VG{C h C 2 ). (98) 
D is computed in 7 operations, and is an approximate solution to all filters of the form: 

R = (A l +A 2 + A 3 +A 4 +A 5 +A 6 +A 7 + A s + c*ONE)»3,whereO <c < 

4 

[123] It can be determined how far off D can be from the correct value of R. The proof 
is omitted for the purpose of brevity, but the answer is R-l <D <#+2. 
[124] After D is computed, the second step is to increase D by 1, so that D is at least as 
big as R. However, one must be cautious on this step: if a byte holds the value of 255 
and it is increased by 1, then most architectures automatically clip the sum of 256 so that 
it becomes 0, which is not what we want. This is remedied by using a "saturated add", 
which will add 1 to each packed byte only if that byte is less than 255. All bytes that are 
255 remain the same. The instruction for this is PADDUSB: 

PADDUSB(D, 1) (99) 
Thus, after 8 instructions, we have computed a value D satisfying R <D <R+3. 
[125] The third step is to determine the correct least significant bits of R. This is done 
by performing the computation of R as it is defined: 

L = CLIP(A } + A 2 + A 3 + A 4 + A 5 + A 6 +A 7 + A 8 + c*ONE)»3. (100) 
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The value c*ONE is a stored constant similar to ONE, so there is no need to perform a 
multiply. Also, most architecture perform the CLIP automatically, so this does not count 
as an extra instruction. In total there are 8 adds and 1 shift to compute L, which holds the 
5 least significant bits of R for each packed byte. 

[126] The fourth step uses L to determine the error term for D. Since we know that D 
is at most 3 more than R, we only need to figure out how much needs to be subtracted 
from D so that it agrees with L in the two least significant bits. This is accomplished by: 

E = CLIP(D -L)& THREE (101) 
As before, the CLIP comes for free on most architectures and THREE is a stored constant 

similar to ONE. In 2 instructions, the error term has been determined. The final step is to 

subtract this error from D to get the final result. This is 1 additional instruction. 

R = D-E (102) 
In total this is 20 instructions, or 19 when c = 0 since an addition can be saved. Pseudo- 
code for the approach using 20 instructions is shown in Table VII, below. 



#define P(a,b) (((a) + (b) + 1 ) » 1 ) 

bl =P(a01,a02); 

b2 = P(a03,a04); 

b3 = P(a05,a06); 

b4 = P(a07,a08); 

cl =P(bl ,b2 ); 

c2 = P(b3,b4); 

d =P(cl ,c2); 

d = ((int)d + 1) < 256 ? d+1 : 255; /* saturated add */ 

e = (aOl + a02 + a03 + a04 + a05 + a06 + a07 + a08 + cc) » 3; 

dd = (d - e) & 0x03; 

xl=d-dd; 

TABLE VII 

[127] The approach of Table VII provides an exact answer in 20 instructions compared 
to 32 of the approach of Table III. In fact, the exact solution provided by Table VII even 
beats the approximate solution of Table VI (22 instructions). However, the approximate 
solution offers significant advantages in many special filter computations. 
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[128] Further refinements can be obtained when certain operands are repeated. For 
instance, consider the following filter: 

(A X +A 2 +2A 2 +2A 5 +2A 1 +4*ONE)»3 
[129] There is no need to do PA VGiA 39 A 3 ) 9 PA VG(A 59 A 5 % or PA VG(A l9 A 7 ) since any 
PA VG of a number with itself is just itself This saves 3 instructions. Further, the 
computation of L can be shortened by first doing (A 3 +A 5 +A 7 ) « 1 and then adding on 
the remaining parts. Hence there are only 5 additions, but 1 additional shift instruction. 
The net savings on L is 2 instructions, and the total net savings is 5 instructions. 
[130] Our results are summarized in Table VIII below. 



Summary of Results for Type 3 FIR Filters. 



Type 3 Filters 


Section 3 Method 


Corrected 
Approximation 
Method 


(A ] +A 2 +... + A s + 4*ONE)»3 


32 


20 


(A x +A 2 + ... + A s + 3*ONE)»3 


35 


20 . 


(A x +A 2 + ... + A 8 + 2*ONE)»3 


34 


20 


(4, +A 2 +... +A s + l*ONE)»3 


35 


20 


(A^+A 2 + ...+A S )»3 


34 


19 


(A i +2A 3 +2A 5 +2A 1 +A 2 +4*ONE)»3 


16 


15 


(A ! +A 2 +A 3 +3A 4 +2A 1 +4*ONE)»3 


19 


17 


(A i +A 2 +A 3 +2A 4 +A 5 +A 6 +A 7 +4*ONE)»3 


27 


19 



TABLE VIII 



B. Corrected Approximation Method Applied to Type 4 Filters 

[131] Following section 4, we define the following packed 64-bit registers, each 

containing 8 data elements of one byte each: 

B X =PA VG(A , A 2 \ B 2 = PA VG(A 3 4 4 ), B 3 = PA VG(A 5 ^ 6 ), B 4 = PA VGiA^z), 
B 5 = PAVG{A 9 A\ 0 ), B 6 =PAVG(A u ^ l2 \ B n = PAVG(A iiy A H ), B s = PAVG(A l5 ^ i6 ), 
C, = PA VG(B X ,B 2 ), C 2 = PA VG(B 3 ,B 4 ), C 3 - PA VG(B 5 ,B 6 ), C 4 = PA VG(B 7 ,B S ), 
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£>! = PA VG(C X ,C 2 ), D 2 = PA FG(C 3 ,C 4 ), D = PA VG(D X J) 2 ) 



(103) 



D is computed in 15 operations, and is an approximate solution to all filters of the form: 

R = (A; + A2 + A3 + A4 +A5 + A6 + A7 + A$ + Ag + Ajo + A]j + A]2 + A/3 + 

A 14 + A/5 + A/o~ + c*ONE)»5 

where 0 <c <8. It can be shown that R-l <D <R+2. 

[132] The remaining steps are more or less the same as done for the type 3 filters: 

PADDUSB(D, 1) 

L = CLIP(A j + A2 + As + A4+A^ + A6 + Ay+As + Ag + Aio + Ajj+A]2 
+ A/3 + A/4 + A/j + A/ 0 ' + c*ONE)»5 

E = CLIP(D -L)& THREE 

R = D-E (104) 

[133] In total this is 36 instructions, or 35 when c = 0 since an addition can be saved. 
The results are summarized in Table IX below. 



Summary of Results for Type 4 FIR Filters. 



Type 4 Filters 


Section 5 
Method 


Corrected 
Approximation 
Method 


{A x + A 2 + ■ . . + A l6 + 8*ONE) » 4 


N/A 


36 ■ 


(A x +AA 2 +6A l +4A A +A s +Z*ONE)»A 


19 


•: 18 


(A ,+/l 2 +2^ 5 +2/l 6 +2v4 7 +2/l 8 +4/l 9 +^3+^ 4 +8*OA'£)» 4 


36 


,24. . 


(A l +2A 2 +2A J +2A A +2A s +2A 6 +2A 1 +2A s +A 9 +i*ONE)»4 


46 




(A l +2A 2 +iA 3 +'iA 4 +3A 5 +2A 6 +A 1 +S*ONE)»4 


36 


:! '■':' 24 y . 



TABLE IX 



[134] Although the invention has been discussed with respect to specific embodiments 
thereof, these embodiments are merely illustrative, and not restrictive, of the invention. 
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For example, although a two-operand SIMD instruction has been primarily discussed, 
techniques and features of the invention may be applicable to other applications where 
the number of operands, or arguments, of a SIMD instruction are not the same as the 
number of variables, or values, in a formula, computation or function to be implemented 
with the SIMD instruction (i.e., a "mismatched" instruction). 
[135] Although the invention has been described with respect to specific SIMD 
instructions to obtain an average of values, any other type of SIMD instruction or 
operation may benefit from the approach of the invention. Although specific operations 
such as addition, subtraction, bitwise AND, bitwise OR, bitwise logical right shift, 
bitwise logical left shift, bitwise exclusive OR, etc., are used in specific embodiments to 
achieve a result, other embodiments may use different operations, or combinations of 
operations, to achieve results. For example, an AND function can be realized by using an 
OR function and complementing, or inverting, the operands and result. Other such 
operational equivalents will be apparent. 

[136] Alternative methods of detecting when the sum of two packed values results in an 
odd number can be employed. Some processors may provide instructions that combine 
multiple operations into compound one or more instructions. Although specific reference 
has been made to a "SIMD" type of instruction, other types of parallel instructions may 
be within the scope of the invention. Although the SIMD instruction has been described 
as a single instruction, other embodiments may use SIMD instructions that occupy more 
than a single instruction's worth of clock cycles, instruction cycles, or the like. 
[137] There are various ways that the invention can be modified from specific 
embodiments described herein to achieve similar results. For example, adjustments to an 
approximate solution are performed as an intermediate step before computing the final 
result so that the approximate solution is no less than the actual solution. One 
modification can be to adjust the approximate solution so that it becomes no larger than 
the actual solution as an intermediate step. Such modifications will be apparent to one of 
skill in the art and are within the scope of the invention. 

[138] Any suitable programming language can be used to implement the routines of the 
present invention including C, C++, Java, assembly language, etc. Different 
programming techniques can be employed such as procedural or object oriented. The 
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routines can execute on a single processing device or multiple processors. Although the 
steps, operations or computations may be presented in a specific order, this order may be 
changed in different embodiments. In some embodiments, multiple steps shown as 
sequential in this specification can be performed at the same time. The sequence of 
operations described herein can be interrupted, suspended, or otherwise controlled by 
another process, such as an operating system, kernel, etc. The routines can operate in an 
operating system environment or as stand-alone routines occupying all, or a substantial 
part, of the system processing. 

[139] Steps can be performed in hardware or software, as desired. Note that steps can 
be added to, taken from or modified from the steps presented in this specification without 
deviating from the scope of the invention. In general, the flowcharts are only used to 
indicate one possible sequence of basic operations to achieve a functional aspect of the 
present invention. 

[140] In the description herein, numerous specific details are provided, such as 
examples of components and/or methods, to provide a thorough understanding of 
embodiments of the present invention. One skilled in the relevant art will recognize, 
however, that an embodiment of the invention can be practiced without one or more of 
the specific details, or with other apparatus, systems, assemblies, methods, components, 
materials, parts, and/or the like. In other instances, well-known structures, materials, or 
operations are not specifically shown or described in detail to avoid obscuring aspects of 
embodiments of the present invention. 

[141] A "computer-readable medium" for purposes of embodiments of the present 
invention may be any medium that can contain, store, communicate, propagate, or 
transport the program for use by or in connection with the instruction execution system, 
apparatus, system or device. The computer readable medium can be, by way of example 
only but not by limitation, an electronic, magnetic, optical, electromagnetic, infrared, or 
semiconductor system, apparatus, system, device, propagation medium, or computer 
memory. 

[142] A "processor" includes any system, mechanism or component that processes data, 
signals or other information. A processor can include a system with a general-purpose 
central processing unit, multiple processing units, dedicated circuitry for achieving 
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functionality, or other systems. Processing need not be limited to a geographic location, 
or have temporal limitations. For example, a processor can perform its functions in "real 
time," "offline," in a "batch mode," etc. Portions of processing can be performed at 
different times and at different locations, by different (or the same) processing systems. 
[143] Reference throughout this specification to "one embodiment", "an embodiment", 
or "a specific embodiment" means that a particular feature, structure, or characteristic 
described in connection with the embodiment is included in at least one embodiment of 
the present invention and not necessarily in all embodiments. Thus, respective 
appearances of the phrases "in one embodiment", "in an embodiment", or "in a specific 
embodiment" in various places throughout this specification are not necessarily referring 
to the same embodiment. Furthermore, the particular features, structures, or 
characteristics of any specific embodiment of the present invention may be combined in 
any suitable manner with one or more other embodiments. It is to be understood that 
other variations and modifications of the embodiments of the present invention described 
and illustrated herein are possible in light of the teachings herein and are to be considered 
as part of the spirit and scope of the present invention. 

[144] Embodiments of the invention may be implemented by using a programmed 
general purpose digital computer, by using application specific integrated circuits, 
programmable logic devices, field programmable gate arrays, optical, chemical, 
biological, quantum or nanoengineered systems, components and mechanisms may be 
used. In general, the functions of the present invention can be achieved by any means as 
is known in the art. Distributed, or networked, systems components and circuits can be 
used. Communication, or transfer, of data may be wired, wireless, or by any other 
means. 

[145] It will also be appreciated that one or more of the elements depicted in the 
drawings/figures can also be implemented in a more separated or integrated manner, or 
even removed or rendered as inoperable in certain cases, as is useful in accordance with a 
particular application. It is also within the spirit and scope of the present invention to 
implement a program or code that can be stored in a machine-readable medium to permit 
a computer to perform any of the methods described above. 
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[146] Additionally, any signal arrows in the drawings/Figures should be considered 
only as exemplary, and not limiting, unless otherwise specifically noted. Furthermore, 
the term "or" as used herein is generally intended to mean "and/or" unless otherwise 
indicated. Combinations of components or steps will also be considered as being noted, 
where terminology is foreseen as rendering the ability to separate or combine is unclear. 
[147] As used in the description herein and throughout the claims that follow, "a", "an", 
and "the" includes plural references unless the context clearly dictates otherwise. Also, 
as used in the description herein and throughout the claims that follow, the meaning of 
"in" includes "in" and "on" unless the context clearly dictates otherwise. 
[148] The foregoing description of illustrated embodiments of the present invention, 
including what is described in the Abstract, is not intended to be exhaustive or to limit the 
invention to the precise forms disclosed herein. While specific embodiments of, and 
examples for, the invention are described herein for illustrative purposes only, various 
equivalent modifications are possible within the spirit and scope of the present invention, 
as those skilled in the relevant art will recognize and appreciate. As indicated, these 
modifications may be made to the present invention in light of the foregoing description 
of illustrated embodiments of the present invention and are to be included within the 
spirit and scope of the present invention. 

[149] Thus, while the present invention has been described herein with reference to 
particular embodiments thereof, a latitude of modification, various changes and 
substitutions are intended in the foregoing disclosures, and it will be appreciated that in 
some instances some features of embodiments of the invention will be employed without 
a corresponding use of other features without departing from the scope and spirit of the 
invention as set forth. Therefore, many modifications may be made to adapt a particular 
situation or material to the essential scope and spirit of the present invention. It is 
intended that the invention not be limited to the particular terms used in following claims 
and/or to the particular embodiment disclosed as the best mode contemplated for carrying 
out this invention, but that the invention will include any and all embodiments and 
equivalents falling within the scope of the appended claims. 
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